Exome Sequencing of Extended Families with Alzheimer's Disease Identifies Novel Genes Implicated in Cell Immunity and Neuronal Function.

OBJECTIVE: Alzheimer's disease (AD) is a neurodegenerative disorder for which more than 20 genetic loci have been implicated to date. However, studies demonstrate not all genetic factors have been identified. Therefore, in this study we seek to identify additional rare variants and novel genes potentially contributing to AD. METHODS: Whole exome sequencing was performed on 23 multi-generational families with an average of eight affected subjects. Exome sequencing was filtered for rare, nonsynonymous and loss-of-function variants. Alterations predicted to have a functional consequence and located within either a previously reported AD gene, a linkage peak (LOD>2), or clustering in the same gene across multiple families, were prioritized. RESULTS: Rare variants were found in known AD risk genes including AKAP9, CD33, CR1, EPHA1, INPP5D, NME8, PSEN1, SORL1, TREM2 and UNC5C. Three families had five variants of interest in linkage regions with LOD>2. Genes with segregating alterations in these peaks include CD163L1 and CLECL1, two genes that have both been implicated in immunity, CTNNA1, which encodes a catenin in the cerebral cortex and MIEF1, a gene that may induce mitochondrial dysfunction and has the potential to damage neurons. Four genes were identified with alterations in more than one family include PLEKHG5, a gene that causes Charcot-Marie-Tooth disease and THBS2, which promotes synaptogenesis. CONCLUSION: Utilizing large families with a heavy burden of disease allowed for the identification of rare variants co-segregating with disease. Variants were identified in both known AD risk genes and in novel genes.

Targeted sequencing of ABCA7 identifies splicing, stop-gain and intronic risk variants for Alzheimer disease.

Several variants in the gene ABCA7 have been identified as potential causal variants for late-onset Alzheimer's disease (LOAD). In order to replicate these findings, and search for novel causal variants, we performed targeted sequencing of this gene in cohorts of non-Hispanic White (NHW) and African-American (AA) LOAD cases and controls. We sequenced the gene ABCA7 in 291 NHW LOAD cases and 103 controls. Variants were prioritized for rare, damaging variants and previously reported variants associated with LOAD, and were follow-up genotyped in 4076 NHW and 1157 AA cases and controls. We confirm three previously associated ABCA7 risk variants and extend two of these associations to other populations, an intronic variant in NHW (P=3.0×10-3) (originally reported in a Belgian population), and a splice variant originally associated in the Icelandic population, which was significantly associated in the NHW cohort (P=1.2×10-6) and nominally associated in the AA cohort (P=0.017). We also identify a 3'-UTR splice variant that segregates in four siblings of one family and is nominally associated with LOAD (P=0.040). Multiple variants in ABCA7 contribute to LOAD risk.

GWAS analysis implicates NF-κB-mediated induction of inflammatory T cells in multiple sclerosis.

To identify genes and biologically relevant pathways associated with risk to develop multiple sclerosis (MS), the Genome-Wide Association Studies noise reduction method (GWAS-NR) was applied to MS genotyping data. Regions of association were defined based on the significance of linkage disequilibrium blocks. Candidate genes were cross-referenced based on a review of current literature, with attention to molecular function and directly interacting proteins. Supplementary annotations and pathway enrichment scores were generated using The Database for Annotation, Visualization and Integrated Discovery. The candidate set of 220 MS susceptibility genes prioritized by GWAS-NR was highly enriched with genes involved in biological pathways related to positive regulation of cell, lymphocyte and leukocyte activation (P=6.1E-15, 1.2E-14 and 5.0E-14, respectively). Novel gene candidates include key regulators of NF-κB signaling and CD4+ T helper type 1 (Th1) and T helper type 17 (Th17) lineages. A large subset of MS candidate genes prioritized by GWAS-NR were found to interact in a tractable pathway regulating the NF-κB-mediated induction and infiltration of pro-inflammatory Th1/Th17 T-cell lineages, and maintenance of immune tolerance by T-regulatory cells. This mechanism provides a biological context that potentially links clinical observations in MS to the underlying genetic landscape that may confer susceptibility.

Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network.

Genetic variation can affect drug response in multiple ways, although it remains unclear how rare genetic variants affect drug response. The electronic Medical Records and Genomics (eMERGE) Network, collaborating with the Pharmacogenomics Research Network, began eMERGE-PGx, a targeted sequencing study to assess genetic variation in 82 pharmacogenes critical for implementation of "precision medicine." The February 2015 eMERGE-PGx data release includes sequence-derived data from ∼5,000 clinical subjects. We present the variant frequency spectrum categorized by variant type, ancestry, and predicted function. We found 95.12% of genes have variants with a scaled Combined Annotation-Dependent Depletion score above 20, and 96.19% of all samples had one or more Clinical Pharmacogenetics Implementation Consortium Level A actionable variants. These data highlight the distribution and scope of genetic variation in relevant pharmacogenes, identifying challenges associated with implementing clinical sequencing for drug treatment at a broader level, underscoring the importance for multifaceted research in the execution of precision medicine.

Genetic overlap between Alzheimer's disease and Parkinson's disease at the MAPT locus.

We investigated the genetic overlap between Alzheimer's disease (AD) and Parkinson's disease (PD). Using summary statistics (P-values) from large recent genome-wide association studies (GWAS) (total n=89 904 individuals), we sought to identify single nucleotide polymorphisms (SNPs) associating with both AD and PD. We found and replicated association of both AD and PD with the A allele of rs393152 within the extended MAPT region on chromosome 17 (meta analysis P-value across five independent AD cohorts=1.65 × 10(-7)). In independent datasets, we found a dose-dependent effect of the A allele of rs393152 on intra-cerebral MAPT transcript levels and volume loss within the entorhinal cortex and hippocampus. Our findings identify the tau-associated MAPT locus as a site of genetic overlap between AD and PD, and extending prior work, we show that the MAPT region increases risk of Alzheimer's neurodegeneration.

Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems.

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the Electronic Medical Records and Genomics Network and the Pharmacogenomics Research Network, has three objectives: (i) to deploy PGRNseq, a next-generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1- to 3-year time frame across several clinical sites; (ii) to integrate well-established clinically validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and to assess process and clinical outcomes of implementation; and (iii) to develop a repository of pharmacogenetic variants of unknown significance linked to a repository of electronic health record-based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site-specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to managing incidental findings, and patient and clinician education methods.

Vascular tone pathway polymorphisms in relation to primary open-angle glaucoma.

AIMS: Vascular perfusion may be impaired in primary open-angle glaucoma (POAG); thus, we evaluated a panel of markers in vascular tone-regulating genes in relation to POAG. METHODS: We used Illumina 660W-Quad array genotype data and pooled P-values from 3108 POAG cases and 3430 controls from the combined National Eye Institute Glaucoma Human Genetics Collaboration consortium and Glaucoma Genes and Environment studies. Using information from previous literature and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, we compiled single-nucleotide polymorphisms (SNPs) in 186 vascular tone-regulating genes. We used the 'Pathway Analysis by Randomization Incorporating Structure' analysis software, which performed 1000 permutations to compare the overall pathway and selected genes with comparable randomly generated pathways and genes in their association with POAG. RESULTS: The vascular tone pathway was not associated with POAG overall or POAG subtypes, defined by the type of visual field loss (early paracentral loss (n=224 cases) or only peripheral loss (n=993 cases)) (permuted P≥0.20). In gene-based analyses, eight were associated with POAG overall at permuted P<0.001: PRKAA1, CAV1, ITPR3, EDNRB, GNB2, DNM2, HFE, and MYL9. Notably, six of these eight (the first six listed) code for factors involved in the endothelial nitric oxide synthase activity, and three of these six (CAV1, ITPR3, and EDNRB) were also associated with early paracentral loss at P<0.001, whereas none of the six genes reached P<0.001 for peripheral loss only. DISCUSSION: Although the assembled vascular tone SNP set was not associated with POAG, genes that code for local factors involved in setting vascular tone were associated with POAG.

Follow-up of loci from the International Genomics of Alzheimer's Disease Project identifies TRIP4 as a novel susceptibility gene.

To follow-up loci discovered by the International Genomics of Alzheimer's Disease Project, we attempted independent replication of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundació ACE data set; 1808 patients and 2564 controls). Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPW1 and FERMT2, respectively. Of these, ZCWPW1 was the only SNP to withstand correction for multiple testing (P=0.000655). Furthermore, we identify TRIP4 (rs74615166) as a novel genome-wide significant locus for Alzheimer's disease risk (odds ratio=1.31; confidence interval 95% (1.19-1.44); P=9.74 × 10(-)(9)).

Genetically meaningful phenotypic subgroups in autism spectrum disorders.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with strong evidence for genetic susceptibility. However, the effect sizes for implicated chromosomal loci are small, hard to replicate and current evidence does not explain the majority of the estimated heritability. Phenotypic heterogeneity could be one phenomenon complicating identification of genetic factors. We used data from the Autism Diagnostic Interview-Revised, Autism Diagnostic Observation Schedule, Vineland Adaptive Behavior Scales, head circumferences, and ages at exams as classifying variables to identify more clinically similar subgroups of individuals with ASD. We identified two distinct subgroups of cases within the Autism Genetic Resource Exchange dataset, primarily defined by the overall severity of evaluated traits. In addition, there was significant familial clustering within subgroups (odds ratio, OR ≈ 1.38-1.42, P < 0.00001), and genotypes were more similar within subgroups compared to the unsubgrouped dataset (Fst = 0.17 ± 0.0.0009). These results suggest that the subgroups recapitulate genetic etiology. Using the same approach in an independent dataset from the Autism Genome Project, we similarly identified two distinct subgroups of cases and confirmed this severity-based dichotomy. We also observed evidence for genetic contributions to subgroups identified in the replication dataset. Our results provide more effective methods of phenotype definition that should increase power to detect genetic factors influencing risk for ASD.

Parkinson disease loci in the mid-western Amish.

Previous evidence has shown that Parkinson disease (PD) has a heritable component, but only a small proportion of the total genetic contribution to PD has been identified. Genetic heterogeneity complicates the verification of proposed PD genes and the identification of new PD susceptibility genes. Our approach to overcome the problem of heterogeneity is to study a population isolate, the mid-western Amish communities of Indiana and Ohio. We performed genome-wide association and linkage analyses on 798 individuals (31 with PD), who are part of a 4,998 member pedigree. Through these analyses, we identified a region on chromosome 5q31.3 that shows evidence of association (p value < 1 × 10(-4)) and linkage (multipoint HLOD = 3.77). We also found further evidence of linkage on chromosomes 6 and 10 (multipoint HLOD 4.02 and 4.35 respectively). These data suggest that locus heterogeneity, even within the Amish, may be more extensive than previously appreciated.

Independent and epistatic effects of variants in VPS10-d receptors on Alzheimer disease risk and processing of the amyloid precursor protein (APP).

Genetic variants in the sortilin-related receptor (SORL1) and the sortilin-related vacuolar protein sorting 10 (VPS10) domain-containing receptor 1 (SORCS1) are associated with increased risk of Alzheimer's disease (AD), declining cognitive function and altered amyloid precursor protein (APP) processing. We explored whether other members of the (VPS10) domain-containing receptor protein family (the sortilin-related VPS10 domain-containing receptors 2 and 3 (SORCS2 and SORCS3) and sortilin (SORT1)) would have similar effects either independently or together. We conducted the analyses in a large Caucasian case control data set (n=11,840 cases, 10,931 controls) to determine the associations between single nucleotide polymorphisms (SNPs) in all the five homologous genes and AD risk. Evidence for interactions between SNPs in the five VPS10 domain receptor family genes was determined in epistatic statistical models. We also compared expression levels of SORCS2, SORCS3 and SORT1 in AD and control brains using microarray gene expression analyses and assessed the effects of these genes on γ-secretase processing of APP. Several SNPs in SORL1, SORCS1, SORCS2 and SORCS3 were associated with AD. In addition, four specific linkage disequilibrium blocks in SORCS1, SORCS2 and SORCS3 showed additive epistatic effects on the risk of AD (P≤0.0006). SORCS3, but not SORCS2 or SORT1, showed reduced expression in AD compared with control brains, but knockdown of all the three genes using short hairpin RNAs in HEK293 cells caused a significant threefold increase in APP processing (from P<0.001 to P<0.05). These findings indicate that in addition to SORL1 and SORCS1, variants in other members of the VPS10 domain receptor family (that is, SORCS1, SORCS2, SORCS3) are associated with AD risk and alter APP processing. More importantly, the results indicate that variants within these genes have epistatic effects on AD risk.

Variations in COL15A1 and COL18A1 influence age of onset of primary open angle glaucoma.

Primary open angle glaucoma (POAG) is a genetically and phenotypically complex disease that is a leading cause of blindness worldwide. Previously we completed a genome-wide scan for early-onset POAG that identified a locus on 9q22 (GLC1J). To identify potential causative variants underlying GLC1J, we used targeted DNA capture followed by high throughput sequencing of individuals from four GLC1J pedigrees, followed by Sanger sequencing to screen candidate variants in additional pedigrees. A mutation likely to cause early-onset glaucoma was not identified, however COL15A1 variants were found in the youngest affected members of 7 of 15 pedigrees with variable disease onset. In addition, the most common COL15A1 variant, R163H, influenced the age of onset in adult POAG cases. RNA in situ hybridization of mouse eyes shows that Col15a1 is expressed in the multiple ocular structures including ciliary body, astrocytes of the optic nerve and cells in the ganglion cell layer. Sanger sequencing of COL18A1, a related multiplexin collagen, identified a rare variant, A1381T, in members of three additional pedigrees with early-onset disease. These results suggest genetic variation in COL15A1 and COL18A1 can modify the age of onset of both early and late onset POAG.

Genome-wide association study of Alzheimer's disease.

In addition to apolipoprotein E (APOE), recent large genome-wide association studies (GWASs) have identified nine other genes/loci (CR1, BIN1, CLU, PICALM, MS4A4/MS4A6E, CD2AP, CD33, EPHA1 and ABCA7) for late-onset Alzheimer's disease (LOAD). However, the genetic effect attributable to known loci is about 50%, indicating that additional risk genes for LOAD remain to be identified. In this study, we have used a new GWAS data set from the University of Pittsburgh (1291 cases and 938 controls) to examine in detail the recently implicated nine new regions with Alzheimer's disease (AD) risk, and also performed a meta-analysis utilizing the top 1% GWAS single-nucleotide polymorphisms (SNPs) with P<0.01 along with four independent data sets (2727 cases and 3336 controls) for these SNPs in an effort to identify new AD loci. The new GWAS data were generated on the Illumina Omni1-Quad chip and imputed at ~2.5 million markers. As expected, several markers in the APOE regions showed genome-wide significant associations in the Pittsburg sample. While we observed nominal significant associations (P<0.05) either within or adjacent to five genes (PICALM, BIN1, ABCA7, MS4A4/MS4A6E and EPHA1), significant signals were observed 69-180 kb outside of the remaining four genes (CD33, CLU, CD2AP and CR1). Meta-analysis on the top 1% SNPs revealed a suggestive novel association in the PPP1R3B gene (top SNP rs3848140 with P = 3.05E-07). The association of this SNP with AD risk was consistent in all five samples with a meta-analysis odds ratio of 2.43. This is a potential candidate gene for AD as this is expressed in the brain and is involved in lipid metabolism. These findings need to be confirmed in additional samples.

A knowledge-driven interaction analysis reveals potential neurodegenerative mechanism of multiple sclerosis susceptibility.

Gene-gene interactions are proposed as an important component of the genetic architecture of complex diseases, and are just beginning to be evaluated in the context of genome-wide association studies (GWAS). In addition to detecting epistasis, a benefit to interaction analysis is that it also increases power to detect weak main effects. We conducted a knowledge-driven interaction analysis of a GWAS of 931 multiple sclerosis (MS) trios to discover gene-gene interactions within established biological contexts. We identify heterogeneous signals, including a gene-gene interaction between CHRM3 (muscarinic cholinergic receptor 3) and MYLK (myosin light-chain kinase) (joint P=0.0002), an interaction between two phospholipase C-ß isoforms, PLCß1 and PLCß4 (joint P=0.0098), and a modest interaction between ACTN1 (actinin alpha 1) and MYH9 (myosin heavy chain 9) (joint P=0.0326), all localized to calcium-signaled cytoskeletal regulation. Furthermore, we discover a main effect (joint P=5.2E-5) previously unidentified by single-locus analysis within another related gene, SCIN (scinderin), a calcium-binding cytoskeleton regulatory protein. This work illustrates that knowledge-driven interaction analysis of GWAS data is a feasible approach to identify new genetic effects. The results of this study are among the first gene-gene interactions and non-immune susceptibility loci for MS. Further, the implicated genes cluster within inter-related biological mechanisms that suggest a neurodegenerative component to MS.

Distribution of COL8A2 and COL8A1 gene variants in Caucasian primary open angle glaucoma patients with thin central corneal thickness.

PURPOSE: One approach to identify genes that contribute to common complex ocular disorders such as primary open angle glaucoma (POAG) is to study the genetic determinates of endophenotypes that are defined by underlying pre-disposing heritable quantitative traits such as central corneal thickness (CCT). Collagen VIII is a major component of Descemet's membrane and studies in mice have indicated that targeted inactivation of the genes encoding the collagen type 8 alpha1 (Col8a1) and collagen type 8 alpha2 (Col8a2) subunits (COL8A1 and COL8A2) results in thinning of the corneal stroma and of Descemet's membrane. The purpose of this study is to evaluate COL8A1 and COL8A2 as candidate genes for thin CCT in human POAG patients. METHODS: 100 Caucasian POAG patients were enrolled in this study. The entire COL8A1 and COL8A2 coding sequence was determined in 8 patients with CCT<513 µm (one standard deviation (36 microns) below the mean (550 microns) and 8 patients with CCT>586 µm (one standard deviation above the mean). Selected COL8A2 exons containing variants of interest were sequenced in the full POAG cohort. Association and quantitative trait analyses were performed. RESULTS: Three patients with CCT less than 513 µm and advanced POAG were found to have missense changes in COL8A2; two patients had a previously identified mutation, R155Q and one had a novel change, P678L (p=0.0035, Fisher's exact test). Missense changes were not found in any of the patients with CCT>513 µm and missense changes in the COL8A1 gene were not found in any patient. One common COL8A2 SNP, rs274754 was also statistically associated with CCT (p=0.018). CONCLUSIONS: In this study we have identified COL8A2 missense changes in a group of Caucasian patients with very thin CCT and advanced POAG. These results suggest that DNA sequence variants in the COL8A2 gene may be associated with thin corneas in some glaucoma patients. Further study of COL8A2 variants in other patient populations, especially those with thinner CCT such as African-Americans would provide further support for a role of COL8A2 in corneal thickness and in glaucoma.

Genetics and pathogenesis of multiple sclerosis.

Multiple sclerosis (MS) is an idiopathic autoimmune neurodegenerative disease. Like many common diseases, MS has a genetic component; however, as with most complex diseases, the genetic architecture may be influenced by heterogeneity, incomplete penetrance, polygenic inheritance, and environmental factors. This clinically complex disease has provided great challenges for geneticists over the years. Although the first consistent genetic association to MS (with HLA-DR*1501) was discovered more than 30 years ago, lack of consistently replicated genetic results has plagued the scientific community. New study design methods (particularly genome-wide associations studies [GWAS]) along with genome project data and larger datasets have allowed several additional MS genes to be identified and consistently replicated. Thus, after many years of frustration, the strong genetic component associated with MS is finally beginning to be characterized.

Examination of tetrahydrobiopterin pathway genes in autism.

Autism is a complex disorder with a high degree of heritability and significant phenotypic and genotypic heterogeneity. Although candidate gene studies and genome-wide screens have failed to identify major causal loci associated with autism, numerous studies have proposed association with several variations in genes in the dopaminergic and serotonergic pathways. Because tetrahydrobiopterin (BH4) is the essential cofactor in the synthesis of these two neurotransmitters, we genotyped 25 SNPs in nine genes of the BH4 pathway in a total of 403 families. Significant nominal association was detected in the gene for 6-pyruvoyl-tetrahydropterin synthase, PTS (chromosome 11), with P = 0.009; this result was not restricted to an affected male-only subset. Multilocus interaction was detected in the BH4 pathway alone, but not across the serotonin, dopamine and BH4 pathways.

Follow-up examination of linkage and association to chromosome 1q43 in multiple sclerosis.

Multiple sclerosis (MS) is a debilitating neuroimmunological and neurodegenerative disease affecting >4,00,000 individuals in the United States. Population and family-based studies have suggested that there is a strong genetic component. Numerous genomic linkage screens have identified regions of interest for MS loci. Our own second-generation genome-wide linkage study identified a handful of non-major histocompatibility complex regions with suggestive linkage. Several of these regions were further examined using single-nucleotide polymorphisms (SNPs) with average spacing between SNPs of approximately 1.0 Mb in a dataset of 173 multiplex families. The results of that study provided further evidence for the involvement of the chromosome 1q43 region. This region is of particular interest given linkage evidence in studies of other autoimmune and inflammatory diseases including rheumatoid arthritis and systemic lupus erythematosus. In this follow-up study, we saturated the region with approximately 700 SNPs (average spacing of 10 kb per SNP) in search of disease-associated variation within this region. We found preliminary evidence to suggest that common variation within the RGS7 locus may be involved in disease susceptibility.

Mutations in the FUS/TLS gene on chromosome 16 cause familial amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal degenerative motor neuron disorder. Ten percent of cases are inherited; most involve unidentified genes. We report here 13 mutations in the fused in sarcoma/translated in liposarcoma (FUS/TLS) gene on chromosome 16 that were specific for familial ALS. The FUS/TLS protein binds to RNA, functions in diverse processes, and is normally located predominantly in the nucleus. In contrast, the mutant forms of FUS/TLS accumulated in the cytoplasm of neurons, a pathology that is similar to that of the gene TAR DNA-binding protein 43 (TDP43), whose mutations also cause ALS. Neuronal cytoplasmic protein aggregation and defective RNA metabolism thus appear to be common pathogenic mechanisms involved in ALS and possibly in other neurodegenerative disorders.